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The  properties  of  laser  sparks  are  investigated  by 
solving  a  simplified,  two-dimensional  energy  equation 
describing  the  steady  state,  subsonic  propagation  of  a  spark 
in  a  channel.  The  propagation  mechanism  is  assumed  to  be 
thermal  conduction  and  thermal  radiation  is  included  as  an 
effective  optically  thin  emission.  The  solution,  obtained 
in  closed  form,  yields  both  axial  and  radial  temperature 
profiles,  as  well  as  the  relationship  between  the  laser 
beam  characteristics,  the  channel  radius  and  the  propagation 
velocity. 

A  solution  of  the  radially  integrated  energy  equation 
is  obtained  by  introducing  the  radial  heat  conduction  param- 
eter A  suggested  by  Raizer.  The  parameter  is  evaluated  by 
means  of  the  solution  to  the  two-dimensional  energy  equation. 
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.  The  model  is  applied  to  air  sparks  of  .15  cm  and 
.50  cm  radius  propagating  at  velocities  up  to  20  cm/sec 
under  the  influence  of  COg  laser  radiation.  The  results 
indicate  a  spark  length  on  the  order  of  1  cm  and  tempera- 
tures near  16,000  °K.  Approximately  half  the  incident 
laser  intensity  is  absorbed  by  the  spark.  Calculation 
of  the  radial  heat  conduction  parameter  verifies  the 
value  of  2.9  used  by  Raizer  for  an  unbounded  spark  of 
,15  cm  radius.  Furthermore,  A  is  found  to  depend  on  the 
spark  radius. 

The  theory  agrees  with  the  experimental  results  for 
stationary  sparks.  However,  it  fails  at  predicting  the 
spark  properties  for  propagation  velocities  of  several 
meters  per  second.  It  is  postulated  that,  at  these 
velocities,  re-absorption  of  thermal  radiation  creates  a 
none qui libri urn  layer  ahead  of  the  spark  front  which  en- 
hances the  absorption  of  laser  radiation  and  increases 
the  propagation  velocity. 
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CHAPTER  I 
INTRODUCTION 

Interest  in  the  interaction  of  matter  with  intense 
laeams  of  coherent  radiation  has  increased  with  the  develop- 
ment of  high  power  lasers.  In  particular,  extensive 
research  is  being  devoted  to  the  investigation  of  dis- 
charges maintained  by  the  absorption  of  laser  radiation. 
These  discharges  are  commonly  referred  to  as  laser  sparks. 

Discharges  are  in  wide  demand  for  technological  and 
diagnostic  applications  and  for  laboratory  studies  of  plasma 
processes.  The  current  methods  for  the  creation  and  main- 
tenance of  discharges  are  based  on  the  release  of  electro- 
magnetic energy  in  a  gas.  Until  recently  discharges,  which 
are  classified  by  the  frequency  of  the  electromagnetic 
field,  ranged  from  the  direct-current  arc  discharge  to  the 
microwave  discharge.  Within  this  frequency  range  the 
transfer  of  electromagnetic  energy  requires  special 
auxiliary  equipment  such  as  electrodes,  conductor  coils 
smd  wave  guides.  As  a  result  the  plasma  is  confined  by 
boundaries  which  can  erode  €ind  introduce  contaminants. 

The  laser  spark  represents  an   extension  of  the  frequency 
range  of  the  generating  field  to  infrared  and  visible 
frequencies.  The  electromagnetic  energy  is  transported  by 


the  laser  beam  thereby  eliminating  the  need  for  special 
equipment.  It  is  possible,  therefore,  to  obtain  pure 
plasmas  burning  in  free  space.  These  characteristics  plus 
the  high  temperatures  attainable  in  laser  sparks  extend 
the  possible  applications  beyond  those  of  lower  frequency 
discharges.  For  example,  a  recent  study  explored  the 
feasibility  of  using  laser  sparks  to  increase  the  enthalpy 
of  the  flowing  gas  in  a  hypersonic  wind  tunnel.  These 
laser-sustained  discharges  also  offer  the  researcher 
interested  in  plasma  physics  a  unique  opportunity  to  study 
an  unconfined  plasma  in  the  laboratory. 

However,  the  occurrence  of  laser  sparks  is  often 
undesirable.  A  laser  beam  of  sufficiently  high  intensity 
can  initiate  a  spark  in  a  gas  or  by  interaction  with  solid 
materials.  The  intensity  of  the  laser  radiation  is 
attenuated  in  passing  through  the  spark  due  to  the  conversion 
of  radiation  energy  to  thermal  energy  which  takes  place  in 
the  spark.  Therefore,  the  intensity  of  laser  radiation 
which  can  be  transported  to  a  target  is  limited  by  the 
threshold  intensity  for  the  creation  of  a  spark.  Exceeding 
this  intensity  results  in  the  creation  of  a  laser  spark 
which  absorbs  a  portion  of  the  incident  radiation  eUiead  of 
the  target. 

Laser  sparks  generally  fall  into  two  categories 
depending  on  whether  the  laser  beam  ignites  the  discharge 
or  merely  maintains  an  externally  initiated  discharge. 
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Ignition  can  occur  in  a  gas  which  is  usually  transparent  to 
laser  radiation  if  the  laser  intensity  creates  sufficient 
ionization  in  a  region  to  initiate  absorption  of  the  laser 
energy. ^"^  Due  to  the  high  intensities  required,  ignition 
has  been  observed  only  with  the  focused  output  of  pulsed 
lasers.  The  initial  breakdown  region  quickly  develops 
into  a  laser  spark.  Shook  waves,  driven  by  the  absorption 
of  laser  radiation,  propagate  outward  from  the  hot  spark 
and  heat  the  surrounding  areas.  In  particular,  the  cold 
gas  ahead  of  the  spark  and  in  the  path  of  the  laser  beam  is 
ionized  and  begins  to  absorb  the  incident  radiation.  Thus, 
a  new  absorbing  layer  is  created  and  the  spark  front 
propagates  in  a  direction  opposite  the  laser  flux  due  to  a 
mechanism  similar  to  that  occuring  in  the  detonation  of 
explosives.  Typical  spark  temperatures  are  10-^  to  10  TC. 
As  the  spark  propagates,  it  is  observed  to  expand  radially 
to  fill  the  cone  formed  by  the  focused  laser  beam. 

In  1969  Bunkin  et  al.^  demonstrated  for  the  first 
time  that  a  laser  intensity  lower  than  that  required  for 
breakdown  could  maintain  an  externally  initiated  discharge. 
They  found  that  Nd-glass  laser  radiation  at  an  intensity 
of  10^  kW/cm^  could  maintain  an  air  spark  which  had  been 
initially  ignited  by  an  electrical  discharge.  This 
Intensity  was  two  orders  of  magnitude  less  than  that 
required  to  initiate  breakdown.  In  this  case  the  incident 
laser  intensity  is  insufficient  for  generating  a  shock  wave. 


1^ 

Thermal  heat  conduction,  aided  by  reabsorptlon  of  thermal 
radiation  emitted  from  the  hot  spark,  is  the  mechanism  by 
which  the  cold  gas  ahead  of  the  spark  Is  Ionized  and 
becomes  opaque  to  the  Incoming  laser  radiation.  The  spark 
propagates  subsonlcally  In  a  direction  opposite  the  laser 
flux.  Typical  spark  temperatures  In  the  subsonic  regime 
are  10,000  -  25,000  °K.^" 

This  mode  of  propagation  suggests  the  existence  of  a 
threshold  laser  Intensity,  such  that  the  absorbed  laser 
energy  Just  compensates  for  the  heat  conduction  and 
radiation  losses,  and  the  spark  remains  stationary.  The 
existence  of  a  threshold  was  demonstrated  In  1970  by 
Generalov  et  al.^'*^  who  used  a  continuous  wave  (cw)  COg 
laser  to  maintain  a  stationary  plasma  In  high  pressure 
argon  and  xenon  gases.  Similar  experiments  have  been 

Q 

conducted  by  Franzen,^  and  recently  Smith  and  Fowler  have 
published  data  on  the  threshold  Intensity  of  COg  laser 
radiation  In  atmospheric-pressure  air. 

The  modeling  of  subsonic  laser  sparks  promotes 
understanding  of  the  creation  and  maintenance  of  these 
discharges.  Unfortunately,  only  a  limited  amount  of 
theoretical  work  has  been  done.  Ralzer,   using  a 
one-dlmenslonal  analysis  and  Introducing  a  thermal 
conduction  parameter  for  radial  heat  losses,  derived 
threshold  laser  Intensities  and  propagation  velocities  for 
sparks  in  atmospheric-pressure  air  under  the  influence 
of  CO9  and  Nd-glass  laser  radiation.  Hall,  Maher  and  Wei, 


5 
resorting  to  analytical  approximations  for  the  thermodynamic 
properties  of  air,  also  used  a  one-dimensional  model  to 
obtain  the  radially  averaged  spark  temperatures  and 
propagation  velocity.  In  both  References  10  and  11,  thermal 
conduction  was  assumed  to  be  the  propagation  mechanism.  A 
recent  one-dimensional  study  by  Jackson  and  Nielsen   has 
focused  on  the  role  of  radiation  in  the  subsonic  propagation 

of  laser  sparks. 

The  purpose  of  this  investigation  is  to  provide  a 
two-dimensional  model  for  the  subsonic  propagation  of  laser 
sparks  by  explicitly  including  the  radial  heat  conduction 
losses.  This  allows  the  calculation  of  both  radial  and 
axial  temperature  profiles  as  opposed  to  temperatures 
averaged  over  the  radius.  Furthermore,  the  solution  is 
obtained  for  a  spark  propagating  in  a  channel  of  arbitrary 
radius  so  that  the  effect  of  boundaries  on  spark  propagation 
can  be  studied.  The  model  is  applied  to  air  sparks  at 
atmospheric  pressure  under  the  influence  of  COg  laser 

radiation. 

The  propagation  mechanism  is  assumed  to  be  thermal 
conduction  with  radiation  serving  only  as  an  energy  loss 
mechanism.  Although  the  flow  velocity  is  restricted  to 
one  dimension,  the  two-dimensional  nature  of  the  temperature 
field  is  retained  in  the  analysis.  By  using  simple  models 
for  the  temperature -dependent  properties  of  the  gas,  a 
closed  form  solution  is  obtained  which  yields  axial  and 


radial  temperature  profiles  as  well  as  the  dependence  of 
the  propagation  velocity  on  channel  radius,  spark  radius 
and  incident  laser  intensity.  The  two-dimensional  model 

also  offers  a  method  for  evaluating  the  radial  heat 

10 
conduction  parameter  introduced  by  Raizer. 
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CHAPTER  II 
THEORETICAL  DEVELOPMENT 

In  this  analysis  the  steady  state,  subsonic 
propagation  of  a  laser  spark  in  a  gas  is  considered. 
The  ignition  process  is  not  considered  since  its  only 
purpose  is  to  provide  free  electrons  to  initiate  inverse 
bremsstrahlung  absorption  of  the  laser  beam.  Therefore, 
the  characteristics  of  the  spark  are  functions  only  of 
the  laser-plasma  interaction  and  not  of  the  ignition 
process. 

The  model  consists  of  a  uniform,  monochromatic  laser 

beam,  of  intensity  S  and  constant  radius  Rj^,  passing  through 

a  channel  of  radius  R^  in  which  there  is  an  absorbing  plasma. 

c 

Energy  is  absorbed  from  the  laser  beam  by  the  plasma  and 
is  transferred  to  the  cold  gas  through  heat  conduction. 
The  channel  wall  is  held  at  a  fixed  temperature  T  =  0  and 
is  assumed  to  absorb  all  radiation  incident  upon  it 
without  re-emission.  The  spark  moves  in  the  channel 
without  distortion  with  a  constant  velocity  u  relative  to 
the  cold  gas,  as  shown  in  Figure  1.  Since  the  flow 
velocities  are  subsonic  the  pressure  gradients  are  small 
and  are  neglected  in  the  analysis.  Furthermore,  the  flow 
kinetic  energy  is  neglected  compared  to  the  thermal  energy. 


8 
since  most  of  the  absorbed  radiation  appears  as  a  temper- 
ature increase  as  opposed  to  directed  fluid  motion.  Viscous 
dissipation  is  also  neglected  since  the  loss  due  to 
dissipation  is  less  than  the  losses  due  to  thermal  radiation 
and  heat  conduction  in  the  region  of  laser-plasma  interaction. 
For  this  model  the  energy  and  continuity  equations  referred 
to  a  co-ordinate  system  fixed  to  the  spark  are 


In  the  above  equations  p     is  the  density,  T  is  the 
temperature,  v  and  v  are  the  axial  and  radial  velocity 
components,  respectively,  c  (T)  is  the  specific  heat 
capacity,   A(T)  is  the  thermal  conductivity,  S  is  the 
intensity  of  the  laser  beam,   Tj'^  (T)  is  the  absorption 
coefficient  at  frequency  )>   and  ^    is  the  net  energy  lost 
by  radiation  per  cm^  per  sec.  The  density  is  assumed  to 
be  related  to  the  temperature  and  pressure  by  the  perfect 
gas  law. 

In  order  to  solve  for  both  the  velocity  field  and  the 
temperature  field,  the  two  momentum  equations  must  be 
included  in  the  analysis.   Clearly,  the  solution  of  these 
four,  nonlinear,  coupled  partial  differential  equations  is 


't^-: 


a  problem  of  considerable  difficulty  unless  simplifying 
assumptions  are  made.  Since  the  temperature  field  is  of 
primary  interest  in  this  analysis,  its  two-dimensional 
nature  is  retained.  On  the  other  hand,  the  velocity 
field  is  simplified  by  assuming  v^  =  0  and  taking  the 
axial  mass  flow  rate  independent  of  radius.  This  one- 
dimensional  approximation  to  the  flow  field  is  often  used 
in  spark  propagation  analysis. -^^"^^  The  equations  are 
then  uncoupled.   In  fact,  Equation  (2.2)  reduces  to 

p  V  =  const.  =  Po  ^  where  j>^    is  the  density  far  ahead 

12 
of  the  spark.  Jackson  and  Nielsen,   on  the  basis  of 

preliminary  results,  predict  little  difference  between 
the  results  obtained  using  the  two-dimensional  flow 
field  and  the  one-dimensional  flow  field  in  radiation- 
dominated  propagation.  For  the  model  considered  here, 
the  error  introduced  by  neglecting  the  radial  deflection 
of  the  flow  should  diminish  as  the  channel  radius  decreases. 

Since  v  must  vanish  at  r  =  0  and  r  =  E^,  the  channel 
r  ^ 

inhibits  the  radial  velocity. 

As  a  result  of  the  one-dimensional  flow  approximation, 
the  energy  equation  reduces  to 


^  i^  =  ^-1?  4-  1  3L  fp  ^U  S?f,    -    <l>  (2.3) 

d    —        


-  : .'  •  ..•  ■■    *.^-.  f.i  t,sf-^-?r-"  -  .*  *pir'K'-V----- 
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where  the  thermal  flux  potential 

e    =     \     A(T)<iT 
o 

has  been  used  to  replace  T  as  the  dependent  variable. 

Equation  (2.3)  can  be  linearized  if  suitable 
approximations  are  made  for  the  temperature  dependence 
of  ?fy,  ^  ,  A  and  c  .  The  first  approximation  is  to 
assume  c  /  A  equal  to  a  constant  which  reduces  cf,    to 
a  constant.  Figure  2  shows  the  variation  of  this  ratio 
for  atmospheric-pressure  air  plotted  from  the  data  of 
Reference  11.  The  definite  temperature  dependance, 
apparent  in  the  figure,  indicates  that  the  assumption  is 
a  shortcoming  in  the  theory.  However,  the  assumption  is 
a  reasonable  first  approximation  in  the  temperature  range 
10,000  -  25,000  °K  which  is  of  greatest  interest  for 
subsonic  laser  sparks.  In  fact,  it  will  be  shown  in 
Chapter  VI  that  this  approximation  is  valid  for  the 
velocities  considered  in  this  analysis.  Although  the 
authors  of  Reference  11  used  analytical  approximations 
for  the  temperature  variations  of  enthalpy  and  thermal 
conductivity,  their  approximations  are  equivalent  to 
assuming  a  constant  c  /  ^  . 

Un-ionized  air  is  transparent  to  COg  laser  radiation, 
which  has  a  wavelength  of  10.6  microns.   However, 
absorption  can  occur  by  inverse  bremsstrahlung  (free-free 
trcmsitions  of  the  electrons)  when  air  becomes  ionized. 


i-,  ,  ^ 
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*  ■ 

In  that  case  the  free  electrons  in  the  presence  of 

scattering  centers  are  accelerated  by  the  electric  field 

of  the  laser  radiation.  The  electron  and  heavy-particle 

temperatures  are  then  equilibrated  by  elastic  collisions. 

Although  both  neutral  particles  and  ions  can  serve 

as  scattering  centers,  absorption  in  the  field  of  neutral 

13 
particles  is  important  only  in  a  very  weakly  ionized  gas. 

At  typical  spark  temperatures  the  ionization  is  appreciable 

so  that  only  absorption  due  to  electron-ion  collisions  is 

considered. 

The  inverse  bremsstrahlung  absorption  coefficient 

for  a  high- temperature,  ionized  gas  is  given  by  the 

13 
modified  Kramers  formula 

A'  =  3.  feSx/0 ^^  (2» 

where  N.  is  the  number  density  of  positive  ions  with 
charge  Z,  N  is  the  number  density  of  electrons,  T^  is 
the  electron  temperature,  ]}     is  the  radiation  frequency 
and  g  is  the  Gaunt  factor.  In  the  single  ionization 
temperature  range,  the  absorption  coefficient  is  approx- 
imately proportional  to  Ng^/T^*.  Therefore,  under  equilibrium 
conditions,  ^f^  will  exhibit  a  rapid  increase  at  temper- 
atures where  the  ionization  becomes  appreciable.  For 
atmospheric-pressure  air  this  temperature  lies  in  the  range 
10,000  °K  to  15,000  °K.    In  fact,  equilibrium  calculations 
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of  ^>;  for  atmospheric-pressure  air  show  a  rapid  Increase 
near  T  =  12,000  °K. 

The  absorption  coefficient  of  air  for  Nd-glass 
laser  radiation,  which  has  a  wavelength  of  1.06  microns, 
consists  of  contributions  from  both  free-free  and  bound- 
free  electron  transitions.  It  also  exhibits  a  rapid 
Increase  at  temperatures  near  12,000  K. 

Consequently,  the  solution  domain  of  the  energy 
equation  is  divided  into  two  regions  separated  by  the 
plane  x  =  0.  For  x<0  the  temperature  is  sufficiently 
low  that  there  is  negligible  absorption  of  laser  radiation, 
i.e.  Xy  =  0.  In  the  region  x  >  0,  the  absorption 
coefficient  is  assumed  to  be  a  constant  so  that 

Sn,  =    So  H(R.-On,  e-'^'^ 

Where  S  is  the  Incident  laser  intensity  and  E{Rj   -  r) 
o  ■^ 

is  the  Heaviside  step  function  defined  to  be  unity  for 

0  -^  r  ^  Rt  and  zero  for  r  >  Rt  .  The  term  "spark  front" 

will  be  used  to  denote  the  plane  x  =  0. 

A  complete  treatment  of  the  emission  and  re-absorption 
of  thermal  radiation  would  require  the  solution  of  a 
complicated  integro-differential  equation.  It  is 
included  in  this  analysis  through  a  simple  phenomenological 
approximation.   In  the  region  x  <  0,  where  the  temperature 
is  comparatively  low  and  decreases  rapidly  with  increasing 
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distance  from  the  spark  front,  the  net  thermal  emission 
is  neglected.  In  the  region  x  >  0,  it  is  included  as  a 
simple  volumetric  loss  term  given  by  the  linear  relation- 
ship 

6  -    rr\  Q 

The  constant  m  is  chosen  to  give  values  of  ^  which 
agree  with  experimentally  determined  radiation  losses 
for  air  at  temperatures  near  15,000  °K. 
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The  analysis  of  Jackson  and.   Nielsen   for  large 
radius  sp»ark8  suggests  that  re-absorption  of  thermal 
radiation,  rather  than  thermal  conduction,  is  the  dominant 
heating  mechanism  for  the  cold  gas  ahead  of  the  spark. 
The  effect  of  neglecting  this  heating,  by  assuming  ^   =  0 
for  X  <  0,  and  the  role  of  radiation  transport  in  spark 
propagation  will  be  discussed  in  Chapter  VI. 

As  a  result  of  the  approximations  discussed  in  the 
preceding  paragraphs,  0,  which  is  the  solution  for  the 
thermal  potential  in  the  region  x  <  0  must  satisfy 

and  ©2  which  is  the  thermal  potential  in  the  region 
X  >  0  must  satisfy 


Ik 


^  ^z      .    ^    r      ^ez\  ^e 


+•-  —  Vr  )-<^  -  rn  0, 


(2.6) 

=    ~  So  H(R,-OK^  e 


-  ^^  X 


The  boundary  conditions  at  x  =  ±  CO      and  at  the  channel 


radius  R^  are 
c 


e,  C-oo,r)=  62(00,0=  e,  Cx,Rc')  =  ei(x,RO=o  (2.7) 

Since  the  temperature  and  the  heat  flux  must  be  continuous 
at  X  =  0,  the  Joining  conditions  are 

0,(0, r)=  0j,CO,r)  (2.8) 


a>^ 


X  =0 


(2.9) 


Equations  (2.7)  -  (2.9) .together  with  the  restriction 
that  0  be  finite,  completely  determine  the  solution. 
There  is,  however,  one  further  condition  which  must  be 
met.  Since  x  =  0  has  been  defined  as  the  plane  where 
appreciable  absorption  begins,  the  temperature  at  x  =  0 
as  determined  from  the  differential  equations  must 
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correspond  to  the  temperature  at  which  the  plasma  becomes 
oi)aque  to  the  laser  radiation.  This  consistency  condition 
is  satisfied  by  requiring  the  thermal  flux  potential, 
averaged  over  the  laser  radius  at  x  =  0,  to  be  equal  to 
0[     ,  the  value  of  the  potential  at  the  ignition 
temperature  Tj^,  i.e. 

^I  =  — ,  \     r  e,  (0,r)dr  (2.10) 

This  equation  results  in  a  relationship  between  the 
propagation  velocity  and  the  required  laser  intensity. 


CHAPTER  III 
SOLUTION  OF  THE  EQUATIONS 

Equations  (2.5)  and  (2.6)  for  the  thermal  potentials 
can  be  solved  in  closed  form  by  the  separation  of  variables 
technique. 

By  assuming 

o 

the  partial  differential  equation  for  ©,  reduces  to  a 
set  of  ordinary  differential  equations 

,2 

2 

where  the  |3^  are  separation  constants. 

The  equation  for  g^^  is  Bessel's  equation,  which  has 
the  solution 

<5n  =  Q-n  To  (p.  r)  +  koYjP.r) 
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J  and  Y  are  the  zero  order  Bessel  functions  of  the  first 
o     o 

and  second  kind,  respectively,  and  the  a^  and  b^  are 
constants.  Since  Yq(0)  =  cP   ,  b^  must  be  zero  for  the 
solution  to  remain  finite  on  the  axis.  The  radial  boundary 
condition  in  Equation  (2.7)  determines  the  separation 
constant.  For  9,  to  vanish  at  r  =  R^,  ^S^  must  equal 
x^^/R„  where  x_„  is  the  n   zero  of  J  . 


on'  c 


on 


The  differential  equation  for  f^^  has  the  solution 


where 


5        -  - 


o^ 
2 


i±J  i  ^^ 


4fil 


and  the  c  and  d^  are  constants.  To  satisfy  the  boundary 
condition  6,  {- oo  ,r)   =  0,  d^^  must  be  zero. 

The  solution  for  ©,  can,  therefore,  be  written  as 


0  =  2-  o-n   To  (>^on  7  )    e 


f  (l-^2fn)X 


(3.1) 


^n=J      i   + 


The  constants  a  are  determined  by  the  joining  conditions. 
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Equation  (2.6)  for  Q^    is  a  nonhomogeneous ,  linear 
partial  differential  equation.  Its  solution  can  be 
written  as 

where  0p  is  a  particular  solution,  and  S\y   is  the 
solution  of  the  associated  homogeneous  equation 


.2 

a  eh 


9^ 


r  3r  V  ^r  J 


(3.2) 


-  m  ev,  =    O 


A  particular  solution  of  Equation  (2.6)  which 
satisfies  the  boundary  conditions  at  r  =  R^  and  x  =  OO  is 


Bp   =    e''^^''  Ik    T.  (xo.^J 


(3.3) 


where  the  b  are  constants.  Substituting  Equation  (3.3) 
into  Equation  (2.6)  yields 


£  -t.  To  (x.o  i  )  =  S„  71^  H  ( R,-  O 


Re 


-t 


fet"  ^*'('<''  +  '')  + 


O.it) 
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Since  the  J«(x^„r/R^)  form  a  complete,  orthogonal 
o  on   c 

set  in  the  range  0  ^  r  ^  R^,  they  can  be  used  to  form 
a  unique  series  expansion  of  any  function  of  r  in  that 
interval.  In  particular,  the  function  S^  "Ky  H(Rj^  -  r) 
can  be  expanded  in  terms  of  an  infinite  series.  Therefore, 
for  Equation  {3 A)   to  be  valid,  the  t^^  must  be  the 
coefficients  of  such  a  series  expansion.  The  coefficients 
can  be  found  by  using  the  orthogonality  property 

0 

where  inm.  is  the  Kronecker  delta  and  J^^  is  the  first 
order  Bessel  fimction  of  the  first  kind.  Multiplying 
Equation  (3>)  by  ^^o^^om^^^o^  *   integrating  from  0  to  R^ 
and  applying  the  orthogonality  condition  gives 


nm. 


■t  = 


R\  T,'(xoO 


The  value  of  the  integral,  obtained  with  the  aid  of 
Equation  (11.3.20)  of  Reference  15,  is  RcVl^^on^L'^^c^'^^on' 
The  resulting  equation  for  the  constants  b^^  of  the 
particular  solution  is 
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1=    — 


>'o^ 


Re  3;'(x„j[  (^Y-  A^  ( ?r^  +  *^  +  m] 


The  solution  of  the  homogeneous  equation  for  Oy^   is 
assumed  to  have  the  form 

Substitution  into  Equation  (3.2)  results  in  the  set  of 
ordinary  differential  equations 


where  the  ^n     are  separation  constants.  The  solutions 
of  these  equations  are  obtained  in  a  manner  analagous 
to  the  procedure  used  for  9,  ,  Imposing  the  boundary 
conditions  F^{  cO  )  =  0   and  G^^^c^  ~  ^  results  in 

©K-  2.c„  T.  (X,  -  )   ••• 


L  = 


-7L{^)'--] 
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where  the  c  are  constants  to  be  determined  by  the  Joining 
conditions. 

Thus,  the  solution  for  62  is  given  by 


^uX 


(3.5) 


tCi-^-))< 


Imposing  the  joining  conditions  given  in  Equations 
(2.8)  and  (2.9)  results  in  a  set  of  algebraic  equations 
for  the  a  and  c  .  Solution  of  these  equations  yields 


the  following  for  O    amd  ©^  : 


e,=  8Ti^5o-r  £  A„  To  (x.„-|  )  e 


t(»-»-Ox 


(3.6) 


e,  =  SKS, 


Re 


-7i;x 


(3.7) 


"ZC^  To(xo.-r  ) 


Re 


e 


iu-^^)y> 


where  the  coefficients  are  given  by 
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^     -  T  (^o<^  Rc^ (3.8) 


(3.9) 


^  ^   c^C^^V^^  ^^^^    ^^  (3.,0) 


The  relationship  between  the  propagation  velocity 
and  the  laser  intensity  is  obtained  by  substituting  S^ 
from  Equation  (3.6)  into  Equation  (2,10).  The  resulting 
equation  is 

ei=  ifc^.  s„  2  —    T,  (x,„  — )         (3.11) 

Since  the  A  are  functions  of  c^  ,  the  propagation 
velocity  is  contained  in  the  summation.  Therefore,  the 
most  convenient  way  to  evaluate  Equation  (3.11)  is  to 
choose  a  value  of  u  and  solve  for  S  .  Equations  (3.6) 
and  (3.7)  can  then  be  evaluated  for  G|  and  0^  . 

Equations  (3.6)  -  (3.11)  are  valid  only  if  each 
infinite  series  converges  uniformly  in  the  interval 
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-  OO  ^  X  £00   and  0  £  r  £  R  .  A  proof  of  the 
uniform  convergence  of  these  series  is  given  in  Appendix  I. 

The  length  of  the  spark  d  can  be  defined  as  the 
co-ordinate  x  where  the  temperature  averaged  over  the 
laser  radius  falls  below  the  ignition  temperature  T^^. 
Although  the  model  does  not  terminate  absorption  of  the 
laser  radiation  at  x  =  d,  the  absorption  of  a  real  spark 
decreases  rapidly  for  x  >  d  due  to  the  rapid  decrease  in 
the  degree  of  ionization.  Therefore,  the  spark  length  is 
useful  in  approximating  the  power  absorbed  by  the  spark. 
Consistent  with  the  exponential  approximation  for  the 
absorption  of  laser  radiation, 

-^  =  I  -  e-"^-"^  (3.12) 

Po 

where  Pa  is  the  power  absorbed  by  the  spark  and  P 
is  the  incident  power. 


CHAPTER  IV 


IMPLIPIED  ONE-DIMENSIONAL  ANALYSIS 


Although  the  solution  in  Equations  (3.6)  -  (3.10) 
is  straightforward,  it  is  cumbersome  due  to  the  summation. 
In  this  chapter  the  equations  are  reduced  to  one  dimension 
by  introducing  a  radial  heat  conduction  parameter.  The 
parameter  is  in  the  form  of  an  infinite  series,  but  the 
spark  properties  are  described  by  simple  equations. 

If  Equations  (2.5)  and  (2.6)  are  averaged  over  the 
laser  radius,  the  resulting  equations  are 


J   e,  de.         A  e, 

—  ex  —  — 


dx^        Jx       Rj 


-   =  0   ,   X  <  O    (4.1) 


2  - 


d  0z  de^      /         A  \    - 

d  x'        d  X      '^       Ry 


(^.2) 
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In  Equations  (^.1)  and  (4.2),  the  radial  heat  conduction 
term  has  been  simplified,  as  suggested  by  Raizer,   by 
setting 

2     /  3e\  A  i 


The  parameter  A  is  assumed  to  be  independent  of  x  and  r 
but  it  is  a  function  of  spark  and  channel  properties. 
For  the  one -dimensional  equations,  the  boundary 
conditions  are 

The  joining  conditions  are  given  by 


d   X  /  X  =0 


The  solutions  of  Equations  (4.1)  and  (4.2)  subject 
to  these  conditions  are 

9^(x^=  G.    e  (4.4) 

—  -  TTv  (  X  -  ci ) 


where 
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^=  \  + 


4A 


.*R^ 


(4.6) 


y  = 


1+  4^.  (mR\  +  A^ 


(4.7) 


(4.8) 


The  relationship  between  the  propagation  velocity  and 
the  incident  laser  intensity  is  given  by 


5o  = 


4^i 


oc  (y  +  2rO 


oc(^'+  l)  ^-  2  7); 


(4.9) 


If  the  losses  due  to  radial  conduction  are  much 

2 
greater  than  the  radiation  losses,  i.e.  A  >>  mR^^  ,  and 

if  the  absorption  coefficient  satisfies  the  relationship 

2  Ap       <<  o<(  2f^+  1),  then  Equation  (4.9)  reduces  to 


5o 


1l 


o('y  C^+  1) 


(4.10) 


10 


Raizer,  ^  neglecting  thermal  radiation  losses  and  the 
attenuation  of  the  incident  laser  beam,  obtained 
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Equation  (4.10)  in  his  one-dimensional  analysis.  For 
large  propagation  velocities,  t      and  }i       both  approach 
unity,  if  the  variation  of  A  with  u  is  neglected.  Then 
Haizer's  result  that  the  laser  intensity  is  proportional 
to  u^  at  large  velocities  is  obtained. 

For  values  of  x  greater  than  the  point  of  maximum 
temperature  for  COp  laser  radiation,  S^     can  be  approx- 
imated by 

0^  ^    Oi    e 


^71^(  X  -  d) 


so  that  d  is  the  length  of  the  spark. 

By  setting  the  derivative  of  O^     in  Equation  (4.5) 
equal  to  zero,  the  location  of  the  maximum  temperature  is 
found  to  be 

oc(^'-0  (l-e-^^'^y 


2  ^. 


Raizer   uses  a  highly  empirical  method  for  calculating 
the  thermal  conduction  parameter.  If  the  axial  variations 
are  neglected  within  a  cylinder  of  radius  R^  with  strongly 
cooled  walls,  the  thermal  potential  will  vary  approximately 
as  J  (xQ,r/RT ) .  When  substituted  into  Equation  (4.3),  this 

2 

distribution  results  in  the  value  A  =  Xq-j^  =5.8.  To 

take  into  account  the  fact  that  the  temperature  is  still 

quite  high  at  r  =  Rj. ,  Raizer  reduces  this  value  by  ^  and 
uses  an  effective  A  =  2.9. 


•  ,. 
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A  more  accurate  determination  of  this  parameter 
can  be  obtained  from  the  two-dimensional  solution 
presented  in  Chapter  III.  By  using  Equation  (3.6)  to 
calculate  9      and  the  radial  derivative,  substituting 
into  Equation  (4.3)  and  evaluating  at  x  =  0  the  following 


\  !.:>''     s'-  ' 


equation  for  A  is  obtained: 

/ \  -  I 

(4.11) 


where  T^^  =  ^Aj^/*qj^)J3^(XqjjRl/Hc)  •  The  convergence  of 
the  infinite  series  is  demonstrated  in  Appendix  I. 


^ 


•*:/V.'.'v' 


CHAPTER  V 
EXAMPLE  CALCULATIONS 

it- 

The  equations  derived  in  Chapters  III  and  IV  have 
been  used  to  study  the  properties  of  sparks  maintained 
in  atmospheric-pressure  air  by  COg  laser  radiation. 

The  variation  of  the  thermal  potential  with  temper- 
ature, obtained  from  the  data  of  Reference  11,  is  shown 
in  Figure  3.  A  value  of  .15  kW/cm  was  used  for  the  value 
of  the  thermal  potential  at  ignition.  This  value  corres- 
ponds to  the  ignition  temperature  Tj^  =  12,000  \   suggested 

10   '  '. 
by  Raizer. 

The  absorption  coefficient  ^^  for  COg  laser  radiation 
in  air  increases  rapidly  for  temperatures  greater  than 
12,000  °K,  due  to  the  rapidly  increasing  degree  of 
ionization. -"-^  It  reaches  a  maximum  value  of  approximately 
.85  cm"^  at  T  «  17,000  °K  and  decreases  to  a  minimum  of 
.38  cm'"^  at  T  =  24,000  °K.  This  behavior  can  be  understood 
by  considering  the  formula  for  7^^,    given  in  Equation  (2.4) 
Beyond  17,000  °K,  first  ionization  is  essentially  complete. 
The  electron  density  approaches  a  plateau  while  the 
temperature  T^  continues  to  increase.  Therefore,   ?fy 
decreases.   For  temperatures  greater  than  24,000  K, 
the  electron  density  again  increases  rapidly  due  to  the 
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.-«^igy;"y- 
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ionization  of  singly  charged  ions.  Consequently,  the 
absorption  coefficient  increases  for  temperatures  greater 
than  24,000  °K.  The  value  K^  =  .7^  cm"^  was  used  in 
the  calculations  for  x  >  0.  It  corresponds  to  an  average 
value  for  the  temperature  range  15,000  °K  to  20,000  °K, 
which  are  typical  values  for  the  temperature  in  the 

si)ark  core. 

The  ratio  c  /A   was  set  equal  to  425  cm  sec/gm. 
From  Figure  1  this  corresponds  to  a  mean  value  in  the 
temperature  range  10,000  °K  to  20,000  \.     The  upstream 
density  was  taken  as  f>^     =  1.3  x  10'-^  gm/cm  . 

The  thermal  emission  parameter  ro  was  assumed  to  be 
175  cm'^.  The  resulting  losses  are  44  kW/cm-^  and 
67  kW/cm^  for  T  =  15,000  °K  and  20,000  °K,  respectively. 
Raizer's  calculations''-^  indicate  losses  which  vary  from 
48  kW/cm^  to  60  kW/cra^  over  the  same  interval,  with  the 
maximum  occurring  at  18,000  °K.  However,  he  uses  only 
half  of  this  emission  in  his  calculations  since  he 
assumes  half  of  the  radiation  is  re-absorbed  in  the  cold 

gas  ahead  of  the  spark  and  not  really  lost.  The  experi- 

1  6  17 
mental  and  theoretical  work  of  Hermann  and  Schade  '   on 

wall-stabilized  nitrogen  arcs  at  atmospheric  pressure 

indicate  radiation  losses  of  50  kW/cm-^  at  20,000  °K. 

Therefore,  the  linear  approximation  with  m  =  175  cm 

adequately  represents  the  losses  in  the  hot  zone  of  the 

air  spark. 


31 
With  these  values  of  the  parameters,  Equations  (3.6)  - 
(3.11)  were  used  to  study  the  propagation  velocities, 
threshold  intensities  and  temperature  profiles  for  air 
sparks  of  radius  .15  and  .50  cm.  Calculations  were 
carried  out  for  propagation  velocities  from  0  to  20  cm/sec. 
The  channel  radius  was  varied  from  R^  =  Rl  to  R^  =  2\> 
A  computer  program  was  written  to  evaluate  the  series. 
It  was  found  that  twenty  terms  were  sufficient  to  insure 
convergence  for  the  range  of  parameters  studied. 

The  variation  of  laser  intensity  with  channel  radius 
for  several  propagation  velocities  is  shown  in  Figure  k. 
For  a  spark  of  .15  cm  radius,  a  decrease  in  channel  radius 
results  in  a  significant  increase  in  the  laser  intensity 
required  to  maintain  the  spark.  The  presence  of  the 
channel  boundary  constricts  the  radial  temperature  profile. 
Therefore,  the  radial  heat  conduction  loss  from  the  spark 

is  increased  as  R^  decreases,  and  a  higher  laser  intensity 

c 

is  required  to  overcome  the  increased  losses. 

The  energy  loss  per  unit  volume  due  to  thermal 
conduction  in  the  radial  direction  is  proportional  to  R^  , 
whereas  the  loss  per  unit  volume  due  to  optically  thin 
radiation  is  independent  of  radius.  Therefore,  thermal 
radiation  becomes  the  dominant  loss  mechanism  for  large 
radius  sparks.  In  that  case  the  presence  of  a  boundary 
should  have  little  effect  on  spark  propagation.  This  is 
demonstrated  in  Figure  '+,  which  shows  that  the  required 
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intensity  for  R,  =  .50  cm  is  nearly  independent  of  channel 

radius . 

Figures  5  and  6  show  the  isotherms,  calculated  from 
Equations  (3.6)  -  (3.11),  for  a  spark  of  .15  cm  radius 
in  a  channel  with  radius  R^  =  1.25  R^.  Figure  5  depicts 
the  threshold  isotherms,  and  Figure  6  shows  the  isotherms 
for  u  =  20  cm/sec.  The  intense  core  and  elongated  struc- 
ture agrees  with  the  qualitative  description  of  laser 

Q 

sparks  given  by  Smith  and  Fowler. 

Figures  7  and  8  show  the  isotherms  for  a  spark  of 
.50  cm  radius  in  a  channel  with  radius  R^  =  1.25  R^.  The 
isotherms  of  Figure  7  correspond  to  the  threshold  case 
while  those  of  Figure  8  correspond  to  a  propagation 
velocity  of  20  cm/sec.  The  blunter  isotherms,  as  compared 
to  the  isotherms  of  Figures  5  and  6,  are  due  to  the 
decreasing  importance  of  radial  heat  conduction. 

The  plots  of  the  isotherms  show  a  rapid  increase  in 
temperature  at  the  spark  front  followed  by  a  slow  decrease. 
The  increase  is  due  to  the  sudden  deposition  of  energy,  as 
a  result  of  the  rapid  increase  in  the  absorption  coef- 
ficient at  T  =  12,000  °K.  The  slow  decline  in  temperature 
for  the  regions  beyond  the  temperature  maximum  is  a 
consequence  of  the  small  absorption  coefficient  for  COg 
laser  radiation.  The  absorption  length  for  GO2  laser 
radiation  in  air  is  Ip    =   TfJ   =  1.^  cm.  Therefore, 
there  is  little  attenuation  of  the  beam  when  it  traverses 
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distances  on  the  order  of  .10  ^^^  .  As  a  result  points 
separated  by  this  distance  absorb  almost  the  same  energy 

from  the  laser  beam. 

Figures  5-8  show  that  the  isotherms  for  propagating 
sparks  are  blunter  at  the  spark  front  when  compared  to 
the  isotherms  for  stationary  sparks.  This  is  a  result  of 
convective  cooling.  In  the  regions  near  x  =  0,  the 
convective  cooling  decreases  with  radius  since  SQ  /  ^^ 
decreases  with  r.  This  Implies  that  the  hotter  regions 
near  the  axis  are  cooled  more  strongly  than  the  outer, 
cooler  regions.  Thus,  convection  tends  to  decrease  the 
radial  variation  of  temperature  near  the  spark  front. 

Representative  values  of  the  physical  properties  of 
laser  sparks,  as  calculated  from  Equations  (3.6)  -  (3.12), 
are  given  in  Table  I.  The  symbol  Qm    denotes  the  maximum 
value  of  the  radially  averaged  potentials,  and  T^^^  is  the 
corresponding  temperature.  For  the  cases  considered  in 
this  analysis  the  spark  is  on  the  order  of  a  centimeter 
in  length  with  average  temperatures  near  16,000  K. 
Approximately  half  the  incident  laser  intensity  is 
absorbed  by  the  spark.  The  higher  intensity  required  to 
maintain  a  propagating  spark,  as  opposed  to  a  stationary 
spark,  results  in  an  increase  in  both  the  temperature 
and  the  spark  length. 

The  variation  of  the  thermal  conduction  parameter  A, 
calculated  for  threshold  conditions,  with  channel  radius 
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is  shown  in  Figure  9.  The  results  indicate  that  A  =  2.9 
for  a  spark  of  .15  cm  radius  in  a  channel  of  .30  cm  radius. 

As  discussed  in  Chapter  IV,  Raizer   obtained  this 
same  value  by  an  approximate  calculation.  The  agreement 
is  only  coincidental  since  his  calculation  indicates  A  is 
independent  of  R^.  In  fact,  the  authors  of  Reference  11 
use  A  =  2.9  for  sparks  with  Rj^  varying  from  .028  cm  to 
1.67  cm.  Figure  9  shows,  however,  that  A  is  a  function 
of  laser  radius.  For  ^q/\   =  2,  the  value  of  A  for  a 
spark  of  .50  cm  radius  is  nearly  twice  the  value  for  a  spark 

of  .15  cm  radius. 

Calculations  of  threshold  laser  intensities  by     ? 
substituting  A  from  Figure  9  into  Equation  (^.9)  gave 
excellent  agreement  with  the  results  of  the  two-dimensional 
analysis.  Although  A  was  found  to  vary  with  propagation 
velocity,  an  error  of  no  more  than  5%   was  incurred  in 
using  the  value  of  A  at  threshold  to  calculate  S^  for 
propagation  velocities  up  to  20  cm/sec. 

The  one-dimensional  solution  can  be  used  to  study 
the  relative  importance  of  each  loss  process  in  subsonic 
spark  propagation.  Dividing  each  loss  term  in  Equation  (^.2) 
by  the  volumetric  energy  absorption  S^H^   exp(-  T^^x) 
results  in  nondimensional  parameters  for  the  losses  due 
to  axial  heat  conduction,  convection,  thermal  radiation    » 
and  radial  heat  conduction.   -  , 
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The  variations  of  these  nondlmensional  parameters 
with  axial  distance  are  shovm  in  Figures  10  and  11  for 
R  =  .15  cm,  R^  =  1.25  Rl  ^'^  propagation  velocities  of 
0  and  20  cm/sec.  These  are  the  same  conditions  used  to 
obtain  the  isotherms  in  Figures  5  and  6.  The  axial 
location  of  the  temperature  maximum  x^^  has  been  used  to 
nondimensionalize  the  axial  distance  from  the  spark  front. 
The  values  of  x^^  are  .132  and  .176  cm  for  u  =  0  and 
20  cm/sec,  respectively. 

Figures  10  and  11  show  that  the  losses  due  to  radial 
heat  conduction  and  to  thermal  radiation  are  approximately 
equal  for  R^  =  .15  cm.  For  x  >  x^^  they  are  the  dominant 
loss  mechanisms.  Axial  heat  conduction  is  important 
only  near  x  =  0,  where  the  rapid  increase  in  temperature 
occurs.  Likewise,  axial  convection  losses,  for  the 
velocities  considered  in  this  analysis,  are  significant 
only  near  the  spark  front,  where   <9e/^>^  is  large. 
Beyond  x  ,  axial  convection  becomes  an  energy  source  due 
to  the  change  in  sign  of  the  axial  gradient  of  the 

temperature . 

Figures  12  and  13  show  the  variations  of  the 
nondlmensional  losses  for  a  spark  of  radius  .50  cm  prop- 
agating at  velocities  of  0  and  20  cm/sec  in  a  channel 
with  R  =  1.25  Rt .  These  are  the  conditions  for  which 

C  L 

the  isotherms  in  Figures  7  and  8  were  plotted.  As  in 
Figures  10  and  11,  the  axial  location  of  the  temperature 
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maximum  has  been  used  to  nondimensionallze  the  axial 
distance.  For  Figures  12  and  13,  x^^^  equals  .129  cm 
and  .231  cm,  respectively. 

Figures  12  and  13  illustrate  the  dominant  role  of 
thermal  radiation  in  the  dissipation  of  energy  for  sparks 
of  large  radius.  Radial  heat  conduction  accounts  for 
only  lOJg  of  the  energy  loss.  Again,  axial  heat  conduction 
and  convection,  for  the  velocities  considered  here,  are 
Important  only  near  the  spark  front. 

It  is  apparent  from  Figures  10  -  13  that  the  laser 
intensity  for  a  propagating  spark  must  exceed  the  threshold 
intensity  so  that  the  additional  loss  due  to  convection 
at  the  spark  front  can  be  overcome.  However,  near  the 
temperature  maximum  convective  losses  are  negligible  so 
that  the  increased  laser  intensity  results  in  a  higher 
maximum  temperatxire  for  a  propagating  spark  than  for  a 
stationary  one.  Beyond  x^^^,  convection  actually  is  a 
heating  term.  The  spark  length  increases  with  propagation 
velocity  due  to  this  additional  source  of  energy  and  due 
to  the  increase  of  laser  intensity  with  propagation 
velocity.  These  results  are  illustrated  in  Figures  5-8 
and  in  Table  I. 


CHAPTER  VI 
DISCUSSION  OF  RESULTS 

Previous  theoretical  studies  of  laser  spark  phenomena 
have  not  probed  the  effect  of  a  channel  boundary  on 
subsonic  spark  propagation.  In  fact,  aside  from  models 

based  on  Raizer's  one-dimensional  approximation,  only  the 

18 
spherically  symmetric  case  has  been  considered.    Con- 
sequently, the  conclusions  drawn  from  this  analysis  cannot 
be  compared  with  previous  theoretical  results. 

The  technology  of  cw  COg  lasers  has  only  recently 
attained  the  power  levels  required  to  maintain  subsonic 
laser  sparks  in  air.  As  a  result,  experimental  investi- 
gations are  scarce,  and  they  have  not  focused  on  spark 
propagation  in  channels.  However,  the  predictions  for 

R  /Rt  «  2  can  be  compared  with  the  results  of  experiments 
C  L 

on  unbounded  sparks  since  the  channel  boundary  is  far 
enough  from  the  spark  that  its  influence  is  negligible. 
A  spark  with  R^/Rl  =  2  will  be  denoted  a  "free"  spark. 

Fowler  et  al..  ^  using  a  focused  COg  laser  beam  in 
atmospheric-pressure  air,  studied  threshold  intensities 
and  propagation  velocities.  They  also  determined  the 
si>atial  variation  of  temperature  and  electron  number 
density  by  a  laser  interferometric  technique.  Their 
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results  indicate  a  threshold  intensity  near  100  kW/cm 
for  a  .15  cm  radius  spark.  This  value  agrees^  well  with 
the  value  of  106  kW/cm^  obtained  from  Figure  4  for  a 
free  spark  of  .15  cm  radius. 

The  isotherms  of  Powler  et  al.  ^  for  an  incident 
power  of  6.2  kW  can  be  compared  with  the  isotherms  of 
Figure  5,  which  correspond  to  an  incident  power  of  7.8  kW. 
The  experimental  results  indicate  an  intense  core,  with  a 
maximum  temperature  of  17,000  °K,  and  an  elongated  struc- 
ture. For  example,  the  maximum  radial  extent  of  the 
1^,000  \   isotherm  is  .10  cm  and  its  length  is  .70  cm. 
From  Figure  5,  the  values  of  the  radial  extent  and  length 
of  the  14,000  °K  isotherm  are  .12  cm  and  .80  cm,  respectively. 
The  authors  of  Reference  19  found  that  thermal  radiation 
was  the  dominant  loss  mechanism  for  Rt  >  .25  cm,  in 

agreement  with  the  results  of  Chapter  V. 

20 
Recently  Keefer  et  al.   studied  the  threshold 

properties  of  laser  sparks  in  air.  They  estimated  the 

threshold  intensity  to  be  120  kW/cm  for  R^  =  .10  cm.  Using 

spectroscopic  techniques,  they  found  isotherms  with  a 

spatial  variation  similar  to  the  variation  reported  by 

Fowler  et  al.  ^  The  maximum  radial  extent  and  length  of  the 

14,000  °K  isotherm  were  found  to  be  .15  cm  and  I.l4  cm, 

respectively. 

In  light  of  the  available  experimental  results,  the 

model  appears  capable  of  predicting  the  threshold  properties 
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of  small  radius  sparks.  It  is  worthwhile  to  consider  the 
simplifying  assumptions  and  the  limitations  they  impose 
on  the  use  of  the  model.      ,•,       ^*";''  r  '   ..j 

The  temperature  dependance  of  c  /  A  was  neglected 
in  the  analysis.  This  approximation  has  no  effect  on 
threshold  calculations  (u  =  0)  since  the  ratio  appears 
only  in  the  convective  loss  term.  For  the  propagation  , 
velocities  considered  in  this  analysis,  convection  is 
significant  only  near  the  spark  front,  as  discussed  in 
Chapter  V.  There  the  temperatures  are  not  much  different 
from  the  ignition  temperature  Tj^.  Therefore,  c^/  ^  can 
be  considered  constant  at  a  value  typical  of  temperatures 
near  12,000  °K. 

The  absorption  coefficient  ?f»;  was  assumed  to  be 
zero  for  x  <  0  and  a  constant  equal  to  .7^  cm"  for 
X  >  0.  The  two  regions  are  separated  by  the  plane  x  =  0 
where  the  temperature  averaged  over  the  laser  reidius 
equals  the  ignition  temperature  T^.  Actually  the  absorp- 
tion coefficient  has  a  maximum  value  of  approximately 
.85  cm"^  at  T  =  17,000  °K  and  decreases  to  .60  cm"  at 
20,000  °K.  Therefore,  temperatures  in  excess  of  17,000  °K 
are  more  difficult  to  achieve  than  the  model  supposes. 

Haizer"^^  has  calculated  maximum  temperatures  for  the 
one -dimensional  case  taking  into  account  the  variation 
of  Tiy,     with  temperature.   His  analysis  shows  that  the 
maximum  temperature  increases  with  propagation  velocity 


and  that  these  temperatures  do  exceed  17,000  TC.  These 
results  support  the  temperature  variation  with  velocity 
indicated  in  Table  I  and  in  Figures  5-8. 

Hall  et  al.^^  compared  the  results  for  a  two-step 
variation  of  ^p  with  a  one-step  model  and  concluded 
that  the  one-step  model  was  a  good  approximation.  Their 
solutions  also  indicate  an  increase  in  temperature 
associated  with  an  increase  in  propagation  velocity. 

Therefore,  the  step- function  formulation  of  ?fy  used 
in  this  analysis  should  not  introduce  appreciable  error 
into  the  calculations  of  S^,  the  required  laser  intensity, 
or  T  ,  the  maximum  of  the  radially  averaged  temperatures. 
However,  an  analysis  including  the  actual  temperature 
variation  of  ^   would  cause  some  suppression  of  temperatures 
greater  than  17,000  \.     For  example,  there  would  be  a 
slight  decrease  in  the  area  within  the  20,000  \   isotherm 
shown  in  Figure  6, 

For  values  of  x  >  d  the  temperature  falls  below 
12,000  °K  and  7i^    decreases  rapidly.  Since  the  model 
retains  the  value  of  .7^  cm"^,  the  temperatures  predicted 
by  the  model  in  that  region  exceed  those  which  actually 
occur.  The  error  introduced  by  this  approximation  is  not 
severe  sl^nce  the  spark  attenuates  the  laser  intensity  to 
approximately  half  its  Initial  value. 

The  analysis  treats  the  transport  of  thermal  radiation 
in  an  approximate  manner.   The  assumption  that  the  radiation 


transport  can  be  represented  by  an  optically  thin  emission 
proportional  to  the  thermal  flux  potential  is  an  oversim- 
plification. Actually,  the  net  energy  loss  due  to  thermal 
radiation  can  only  be  found  by  solving  the  radiative  transfer 
and  energy  equations  simultaneously.  As  discussed  in 
Chapter  V,  the  role  of  radiation  increases  with  spark  radius. 
As  a  consequence,  sparks  of  large  radius  may  require  a  more 
carefiQ  treatment  of  the  radiation  transport. 

The  analysis  presented  here  has  been  limited  to  prop- 
agation velocities  of  20  cm/sec  or  less.  However, 
experiments^*® '^^  indicate  that  the  combustion  front 
velocities  can  exceed  several  meters  per  second.  The  prop- 
agation of  sparks  at  such  high  velocities  by  thermal 
conduction  alone  leads  to  excessively  high  temperatures. 
For  example,  the  solution  of  Equations  (3.6)  -  (3.11)  for 
a  free  spark  of  radius  .50  cm  propagating  at  1  ra/sec 
yields  a  required  laser  intensity  of  7^0  kW/cm  .  This 
intensity  results  in  a  value  of  T^^^  exceeding  40,000  °K, 
which  is  a  factor  of  two  greater  than  experiments 
indicate.^*®'-"-^ 

Due  to  the  thermal  expansion  of  the  spark,  the 
velocity  in  the  laboratory  frame  can  exceed  u  which  is  the 
velocity  measured  relative  to  the  cold  gas  ahead  of  the 
spark. ^»-*-°»-'-®  The  magnitude  of  the  difference  can  only  be 
obtained  by  solving  the  full  two-dimensional  flow  equations. 
It  is  doubtful,  however,  that  thermal  expansion  could 
account  for  the  large  discrepancy. 
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The  discrepancy  between  observed  and  predicted 
velocities  has  led  to  a  closer  examination  of  the  role 
of  thermal  radiation  in  spark  propagation.  Several 
researchers^^-^^'^^  have  noted  that  radiation  emitted  by 
the  hot  spark  can  be  re-absorbed  by  the  cold  gas  ahead  of 
the  spark.  Jackson  and  Nielsen^^  have  shown  that,  for 
large  radius  sparks,  heating  of  the  cold  gas  by  re- 
absorption  can  far  exceed  the  heating  due  to  thermal 
conduction.  Their  results  predict  propagation  velocities 
of  several  meters  per  second  with  temperatures  ^^   of 
approximately  20,000  °K.  However,  the  velocities  they 
obtain  are  still  a  factor  of  five  less  than  the  experi- 
mental  results  they  use  for  comparison. 

The  experiment  of  Keefer  et_al.^°  suggests  that  highly 
nonequilibrium  processes  may  be  occurring  in  the  cooler 
plasma  regions.  They  found  that  the  dominant  radiation 
in  the  sheath  surrounding  the  hot  plasma  core  came  from 
the  1^*  negative  system  of  ^^ ,     Furthermore,  radiation 
from  the  2^^  positive  system  of  Ng,  usually  present  in 
high-temperature  air,  was  conspicuously  absent.  The 
anomalous  Ng"^  radiation  indicates  that  nonequilibrium 
processes,  such  as  phot oionizat ion,  may  be  occurring. 
This  result  suggests  another  explanation  for  the 
high  propagation  velocities.  The  thermal  radiation 
escaping  the  plasma  core,  rather  than  raising  the  temper- 
ature of  the  cold  gas  as  a  whole  as  assumed  in  Reference  12, 


^3 
could  lead  instead,  through  photolonization  and  inverse 
bremsstrahlung  absorption,  to  a  nonequilibrium  region 
ahead  of  the  spark.  A  discussion  of  this  phenomenon 
and  its  effects  on  spark  propagation  is  given  in  Appendix  II 
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CHAPTER  VII 
CONCLUSIONS 

The  two-dimensional ^energy  equation  describing  spark 
propagation  in  a  channel  has  been  used  to  investigate  the 
properties  of  subsonic  laser  sparks.  The  closed-form 
solution  yields  both  axial  and  radial  temperature  profiles, 
as  well  as  the  dependence  of  the  propagation  velocity  on 
the  incident  laser  intensity,  laser  radius  and  channel 

radius .  . 

The  model  was  applied  to  .15  cm  and  .50  cm  radius 
sparks  driven  by  COg  laser  radiation  at  velocities  up  to 
20  cm/sec  in  atmospheric-pressure  air.  The  results 
indicate  the  channel  has  a  significant  influence  on  spark 
propagation  for  the  .15  cm  spark  when  the  channel  radius 
is  less  than  twice  the  spark  radius.  For  sparks  of  .50  cm 
radius,  the  effect  of  the  channel  is  negligible. 

For  the  cases  considered,  the  spark  length  ranges 

from  .3^   cm  to  1.26  cm.  The  maximum  value  of  the  radially 

averaged  temperature  varies  from  13,000  °K  to  17,000  ^. 

The  ratio  of  the  power  absorbed  by  the  spark  to  the 

incident  power  ranges  from  22^  to  61^.  For  a  fixed  R^ 

and  H  ,  these  quantities  increase  with  increasing  laser 
c 

intensity. 
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For  the  .15  cm  spark,  radial  heat  conduction  and 
thermal  radiation  losses  are  approximately  equal.  However, 
for  the  .50  cm  spark,  radiation  is  the  dominant  loss 
mechanism.  In  both  cases,  axial  convection  and  axial  heat 
conduction  are  significant  only  near  the  spark  front. 

The  predicted  threshold  intensity  and  isotherm 
variations  for  the. 15  cm  spark  agree  with  the  available 
experimental  data.  However,  the  model  fails  to  predict 
the  high  propagation  velocities  experimentally  observed. 
In  order  to  maintain  such  high  velocities,  thermal  radi- 
ation must  play  a  significant  role  in  the  heating  of  the 
cold  gas  ahead  of  the  spark.   ? 

An  equation  for  the  thermal  conduction  parameter  A, 
introduced  by  Raizer,  has  been  obtained  from  the  two- 
dimensional  solution.  The  evaluation  of  A  for  a  free  spark 
of  .15  cm  reidius  confirms  the  value  of  2.9  used  by  Raizer, 
Furthermore,  A  was  found  to  depend  on  channel  and  laser 
radius . 

In  summary,  the  two-dimensional  model  presented  in 
this  analysis  can  be  used  as  a  reasonable  approximation 
for  the  properties  of  laser  sparks  at  threshold  and  at 

low  propagation  velocities.  In  those  cases,  the  assumption 
that  thermal  radiation  serves  only  as  a  loss  mechanism 
is  adequate.  For  high  propagation  velocities,  a  more 
careful  treatment  of  the  radiation  transfer  is  required. 


In  spite  of  the  rapidly  expanding  theoretical  and 
experimental  efforts,  understanding  of  subsonic  spark 
phenomena  is  far  from  complete.  In  particular,  the 
.observed  propagation  velocities  have  eluded  a  satisfactory 
explanation.  Further  analysis  in  two  areas  would  provide 
greater  insight  into  the  physical  processes  occurring 
in  subsonic  sparks.  A  detailed  study  of  the  radiation 
transport,  including  nonequilibrium  effects,  would  provide 
a  more  complete  knowledge  of  the  propagation  mechanism. 
And  a  study  of  the  two-dimensional  flow  field  around  the 
spark  would  lead  to  a  more  accurate  determination  of  the 
role  convection  plays  in  the  energy  balance. 
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Figure  2.  Variation  of  c  /A   with  Temperature 
for  Atmospheric-pressure  Air. 


4-9 


CO 


<D 


o 

-p 
o 
Pi 


CD 


^ 

^ 


E 
u 


« 

St 
bO 


50 


Figure  4.  Incident  Laser  Intensity  Required  to  Maintain 
a  Spark  as  a  Function  of  Channel  Radius. 
Solid  curves  are  for  Rt  =  .15  cm  and  broken 
curves  are  for  R-  =  .50  cm.  Curves  1,  2  and  3 
correspond  to  velocities  of  0,  10  and 
20  cm/sec,  respectively. 
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Figure  7.   Isotherms  for  a  Spark  with  R^  =  .50  cm, 
R 


=  1.25  R^  and  u  =  0. 


-C    L -■   (1)  '^  =  0' 

(2)  T  =  8,000  °K,  (3)  T  =  10,000  °K, 

(J^)  T  =  13,000  °K. 
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Figure  8.  Isotherms  for  a  Spark  with  R^^  =  .50  cm, 

Rq  =  1.25  \   and  u  =  20  cm/sec.   (1)  T  =  0, 
(2)  T  =  10,000  °K,  (3)  T  =  14,000  °K, 
(4)  T  =  17,000  °K. 


f.    •  .::■ 


55 


12 


6  - 


, 

, 

\ 

A 

»                       1 

1                       1 

x.o 


1.5 


2.0 


Figure  9.  Variation  of  the  Thermal  Conduction  Parameter 
with  Rq/Rl.  Solid  curve  is  for  R-  =  .15  cm 
and  broken  curve  is  for  Rr  =  .50  cm. 
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Figure  10. 


Axial  Variation  of  the  Volumetric  Energy- 
Loss  Terms  for  B^  =  .15  cm,  R^  =  1.25  Rj^ 
and  u  =  0.   (1)  Radial  heat  conduction, 
(2)  Thermal  radiation,  (3)  Axial  heat 
conduction. 
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Figure  11. 


Axial  Variation  of  the  Volumetric  Energy  Loss 
Terms  for  R^  =  .15  cm,  R  =  1.25  Rt  and 
u  =  20  cm/sec.   (1)  Radial  heat  conduction, 
(2)  Thermal  radiation,  (3)  Axial  heat 
conduction,  (^)  Convection. 
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Figure  12.  Axial  Variation  of  the  Voliunetric  Energy- 
Loss  Terms  for  R^  =  .50  cm,  R^  =  1.25  R^ 
and  u  =  0.   (1)  Thermal  radiation, 
(2)  Radial  heat  conduction,  (3)  Axial 
heat  conduction. 
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Figure  13.  Axial  Variation  of  the  Volumetric  Energy- 


Loss  Terms  for  Rt  =  .50  cm,  R  = 

Li  C 


1.25  R, 


and  u  =  20  cm/sec.  (1)  Thennal  radiation, 
(2)  Radial  heat  conduction,  (3)  Axial  heat 
conduction,  (k)   Convection. 
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Table  I 

Physical  Properties  of  Laser  Sparks 
of  .15  cm  and  .50  cm  Radius 


\ 

u 

^m 

^m 

d 

Pa/^o 

cm 

cm/sec 

kW/cm 

10^  \ 

cm 

Sl  = 

.15  cm 

.15 

0. 

,2k 

15. 

.83 

.46 

.15 

10. 

.28 

16. 

1.02 

.53 

•  .15 

20. 

.32 

17. 

1.24 

.60 

.30 

0. 

.21 

1^. 

.66 

.39 

.30 

10. 

.25 

15. 

.94 

.50 

.30 

20. 

.31 

17. 

1.26 

.61 

■  ■  • 

«L  = 

.50  cm 

.50 

0. 

.19 

13.5 

.50 

.31 

.50 

10. 

.23 

1^.5 

.84 

.46 

.50 

20. 

.30 

16.5 

1.26 

.61 

1.00 

0. 

.17 

13. 

.34 

.22 

1.00 

10. 

.22 

14. 

.74 

.42 

1.00 

20. 

.29 

16. 

1.21 

.59 

APPENDIX  I 
CONVERGENCE  OF  THE  SERIES  SOLUTION 

The  two-dimensional  energy  equation  describing 
subsonic  spark  propagation  was  solved  in  Chapter  III 
in  terms  of  infinite  series.  For  the  solution  to  be 
valid,  the  series  must  be  shown  to  converge  uniformly 
for  every  point  (x,r)  in  the  solution  domain  ~   QO    ^     x  £  cO  , 

0  £  r  £  Rg. 

The  following  properties  of  Bessel  functions  and 
their  zeroes  will  be  used  in  proving  the  convergence  of 

21 
the  series: 

I  Jq(z)  111  <A1.1) 

I  Ji(z)  1^1  (A1.2) 

J^(Xq^)  ^     0  (A1.3) 

lim  Xqjj  «  (n  -  i)  ^  (Al» 

lim  J,(z)  »  .798  z"*  cos(z  -  3  ^  A)   (A1.5) 

lim  x^^  Ji2(x^^)  =  .637  (A1.6) 

The  uniform  convergence  of  a  series  of  function  can 

22 
be  proved  by  the  Weierstrass  M-Test.    It  states  that  the 

series   ^  f  (x,r)  converges  uniformly  in  an  interval  of 

the  variables  x  and  r  if  |fn(x,r)  |  £   M^  in  that  interval 

and   21  M  converges. 
n.   " 
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In  both  Equations  (3.6)  and  (3.7),  the  exponential  term 
and  Bessel  function  which  multiply  the  A^,  B^^  and  C^^  are 
bounded,  in  absolute  value,  by  one.  For  the  series  in 
Equation  (3.6),  M^  can  be  taken  as 


^.= 


6l'  j^;  i^.-hu)  X,.  j;'(x,j 


-  \ 


since  \  A  1  1  M  .  From  the  definitions  of  t^    and  t 
and  Equation  (A1.4),  .  „, 

/       ZfTfl 
ft— *oo  l\-»«»  OtHc- 

Therefore, 

lim  M^  =  Dn"^ 
n  -►  oo  ■  i- 

where  D  is  a  positive,  finite  constant  for  R^  <  <?o  . 
Since  the  series  'Z  n"^  converges  and  since 

rv 

lim  M^/n"^  =  D  <  oo 

22 
then   £.  M  converges  by  comparison.    Consequently,  the 

infinite  series  in  Equation  (3.6)  converges  uniformly.  The 

convergence  of  Equation  (3.11)  follows  by  comparison  with 

the  series    2  A  . 
n.    n 

At  first  glance,  the  two  series  in  Equation  (3-7) 
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appear  to  have  singularities  at  values  of  velocity  such  that 

However,  when  this  condition  is  satisfied,  the  exponential 

factors  for  the  two  series  are  equal.  The  two  series  can 

then  be  written  as  one  series  with  coefficients 

D  =  B  +  C  .  When  this  is  done,  the  singularity  disappears, 
n    n    n 

By  the  Weierstrass  M-Test,  the  series  in  Equation  (3.7) 
are  uniformly  convergent  if  the  series   21  \^^\  ^^^ 

2.  |c  I  converge.  The  convergence  is  easily  shown  by 
comparison.  Since 

rx  -♦oo 


""    I  =n  I  /  h"  I  " 


ft  -♦oo 


and   Z  |a^  I  converges,  then   2.  Pn  ^^    2  ^n 
converge . 

The  infinite  series   2.  ^^^^^  T^  in  Equation  (^.11) 
can  be  shown  to  converge  by  comparison  with  the  series 
Z_  n'^^^.     By  applying  Equations  (A1.4)  -  (A1.6)  it  can 
be  shown  that 


Xo\  T.   _   P-c  cos     [x,^    Re  -  ^  / 


a-»tfo      a"^^^     8lrR».  ^-*oo  -tX  *'^ 

=  0 

Therefore,  the  series   21  ^on^  "^n  converges. 


••■"■  :  ■  -IT' 


APPENDIX  II 
-  .    EFFECT  OF  RADIATION  ON  SUBSONIC  SPARKS  ,^^  v  ;  ^ 

It  was  determined  in  Chapter  VI  that  spark  propagation 
at  velocities  of  several  meters  per  second  by  thermal  con- 
duction results  in  temperatures  which  exceed  the  experimentally 
observed  temperatures.  Therefore,  it  has  been  postulated 
that  re-absorption  of  thermal  radiation  plays  a  dominant 
role  in  the  propagation  of  laser  sparks  at  high  velocities. 

Jackson  and  Nielsen^^  obtained  propagation  velocities 
of  several  meters  per  second  by  including  radiative  trans- 
fer in  a  one-dimensional  model.  Their  solution  was  based 

on  calculations  of  the  radiative  emission  and  absorption 

20 
for  equilibrium  air.  The  recent  results  of  Keefer  et_al. 

suggest,  however,  that  highly  nonequilibrium  processes  are 
occurring  in  the  cooler  plasma  regions.      ^     ; 

Except  for  two  limiting  cases,  the  solution  of  the 
energy  equation  including  radiative  transfer  poses  a  problem 
of  considerable  difficulty.  In  the  optically  thick  limit, 
where  the  mean  free  path  of  the  thermal  radiation  ly  is 
much  less  than  some  characteristic  dimension  of  the  spark  L, 

the  radiative  energy  flux  can  be  written  as   ^^  ^  ^'  *'^®^® 

23 
\^   (T)  is  the  radiative  analogue  of  the  thermal  conductivity. 

In  that  case,  the  effect  of  radiative  transport  is  merely  an 


6it 
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increase  in  the  thermal  conduction  coefficient.  In  the 
optically  thin  limit,  where  iy»  L,  the  plasma  acts  as 

a  volume  radiator,  and  the  rate  of  energy  loss  per  unit 

23 
volume  can  be  written  as  a  function  of  temperature.    Both 

limiting  cases  are  valid  only  where  local  thermodynamic 
equilibrium  exists. 

The  model  described  in  Chapter  II  assumes  that  radia- 
tion can  be  treated  as  an  effective,  optically  thin  emission. 
This  assumption  is  adequate  for  small  radius  sparks  propaga- 
ting at  low  velocities.  In  that  case,  energy  transport  by 
thermal  radiation  is  no  greater  than  energy  transport  by 
thermal  conduction. 

For  laser  sparks  in  which  radiation  dominates  the 
energy  transport,  the  radiative  transfer,  unfortunately,  is 
neither  optically  thick  nor  optically  thin.  The  radiation 
escaping  the  hot  plasma  core  lies  in  the  vacuum  ultra- 
violet (vuv)."^^  In  the  core  region,  where  there  is  a  high 
degree  of  ionization,  the  plasma  is  optically  thin  to  vuv 
radiation.  However,  in  the  cooler  regions  this  radiation  is 
strongly  absorbed  by  the  neutral  species. 

Jackson  and  Nielsen^^^  foxind  that  most  of  the  emission 
from  the  plasma  core  lies  in  the  energy  range  from  1^.5  to 
20  ev,  due  mainly  to  electron  recombination  with  N  and 
o"*"  ions.  Photons  within  that  range  possess  energy  exceeding 
the  ionization  potentials  E.  of  the  neutral  species  of  air, 
E^   being  15.6  ev,  12.1  ev,  14.5  ev  and  13.6  ev  for  Ng,  Og, 
N  and  0,  respectively.    The  phot oionizat ion  cross  sections 
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for  photons  whose  energy  exceeds  E^^  are  high.  For  Ng  and 
Op  these  cross  sections  are  on  the  order  of  20  Mb,  where  a 
megabam  (Mb)  equals  lO'-*-^  cm^.  For  N  and  0  the  cross  sec- 
tions  are  approximately  8  Mb. 

The  absorption  coefficient  for  the  photons  equals  the 
sum,  over  the  absorbing  species  present,  of  the  product  of 
the  cross  section  for  photoionization  and  the  number  density. 
Equilibrium,  atmospheric-pressure  air  at  12,000  °K  is  com- 
posed  primarily  of  neutral  nitrogen  and  oxygen  atoms. 
The  combined  neutral  atom  densities  are  on  the  order  of 
5  X  10^^  cm"^.  A  cross  section  of  8  Mb  results  in  a  absorp- 
tion coefficient  for  the  photoionizing  radiation  of  k   cm" 
and  a  photon  mean  free  path  of  Ly    =  .25  cm.  Consequently, 
the  radiation  is  strongly  absorbed  in  air  at  temperatures 
near  12,000  °K. 

This  strong  absorption  indicates  the  radiation  transport 
in  the  region  x  <  0  might  be  adequately  represented  by 
the  diffusion  approximation   Ar^  '^ -     Jackson  and  Nielsen 
found,  however,  that  this  approach  yielded  values  of  the 
energy  deposition  by  re-absorption  which  were  a  factor  of 
forty  smaller  than  their  numerical  calculations.  The  large 
discrepancy  is  due  to  the  rapid  temperature  variation,  near 
the  spark  front,  over  distances  on  the  order  of  i^f   . 
As  a  result  the  radiation  field  is  highly  anisotropic  and 
the  diffusion  approximation  is  invalid. 

The  effect  of  the  absorption  at  temperatures  near 
12,000  °K  is,  initially,  an  increase  in  the  concentrations 
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of  N"*",  O"*"  and  electrons  above  their  equilibrium  values.  If 
the  relaxation  time  is  short  enough,  the  concentrations  will 
equilibrate  to  values  characteristic  of  some  higher  tempera- 
ture. The  net  effect  of  the  radiation  is  then  an  equilib- 
rium heating  of  the  cold  gas,  as  assumed  by  Jackson  and 
Nielsen. "^^  One  reaction  which  is  essential  for  equilibration 
is  the  electron-recombination  reaction 

N"*"  +  e  +  M  ►  N  +  M  (A2.1) 

where  M  is  an  atomic  or  molecular  species.  At  electron 
temperatures  near  12,000  °K,  the  reaction  rate  constants 
for  such  reactions  are  on  the  order  of  lO"  ^  cm  /sec. 
There  is,  however,  the  competing  reaction    ,. ,  ■' ■ 

N"^  +  N  +  M  — ►  Ng"*"  +  M  ..  .  ^..  .   .  (A2.2) 

for  the  excess  n"*".  The  constant  for  this  reaction  has 
not  been  accurately  determined.  An  estimate  for  tempera- 
tures near  12,000  °K  is  10'^^  -  ^  cm^/sec.^"^  For 
equilibrium,  atmospheric-pressure  air  at  12,000  K  the 
nitrogen  atom  number  density  is  approximately  six  times 
the  electron  number  density, "'■  so  that  reaction  (A2.2)  may 
proceed  at  a  much  faster  rate  than  reaction  (A2.1).  Conse- 
quently, a  nonequilibrium  concentration  of  Ng  may  exist 
outside  the  hot  plasma  core.  In  fact,  investigators  ' 
have  reported  strong  N^^  radiation  outside  the  core  region. 
The  role  of  reaction  (A2.1)  in  the  relaxation  scheme  increases 
with  increasing  electron  number  density.  The  relative 
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importance  of  reactions  (A2.1)  and'(A2.2)  may  change, 
therefore,  as  the  departure  from  equilibriiwa  increases. 

At  6,000  °K  equilibrium  air  consists  mostly  of  N2, 

1^ 
N  and  0,  in  approximately  equal  concentrations. 

Therefore,  if  the  distance  in  which  the  temperature  drops 

from  12,000  °K  to  6,000  °K  is  on  the  order  of  it,    =  .25  cm 

or  less,  some  radiation  can  penetrate  into  the  cooler 

regions  and  directly  photolonize  Ng  by  the  reaction 

« 

N^  +  hP   — ►  Ng"*"  +  e,   h>^  >  15.6  ev       (A2.3) 

It  is  interesting  to  note  that  the  molecular  ions 
produced  by  reaction  (A2.3)  f or  h  v*  >  18.7  ev  are  almost 
equally  distributed  among  the  ground  state  and  the  first 
two  excited  states,  A  ^  fj^  and  B  ^  1.  ^ }       It  is  the 
transition  from  B  ^  21u'^  ^°  ^^®  ground  state  X   Xg 
which  produces  the  1^^  negative  band  of  Ng"*"  observed 
experimentally  in  References  8  and  20. 

A  careful  study  of  the  reactions  which  occur  sub- 
sequent to  photoionization  is  necessary  for  determining 
the  conditions  at  the  spark  front.  Reactions  related  to 
electron  loss  or  gain  are  especially  important  since  the 
absorption  of  laser  radiation  is  directly  related  to  the 
electron  number  density.  A  major  electron  loss  process 
is  the  dissociative  recombination  of  Ng 

N^"^  +  e  ►  2N  (A2.^) 
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Its  rate  constant  is  on  the  order  of  10"'  cm-^/sec. 
Three-body  recombination 

X"^  +  e  +  M  — ^  X  +  M  (A2.5) 

where  X  is  N,  0  or  Ng,  also  plays  a  significant  role. 
The  reaction  rate  constant  for  (A2.5)  at  electron 
temperatures  near  6,000  °K  is  typically  5  x  10'  ^  cm  /sec 
when  M  is  an  atomic  or  molecular  species  and 
5  X  10"^^  cm^/sec  when  M  is  an  electron. ^"^  Electron 
attatchment  to  oxygen  atoms  and  molecules  may  also  be 

important  since  the  negative  ions  can  react  with  other 

29  30 
species  to  form  stable  compoimds.  ^»^ 

If  the  recombination  reactions  are  fast  enough,  the 

electron  number  density  will  relax  to  its  equilibrium  value 

before  the  electrons  can  be  appreciably  heated  by  the  laser 

radiation.  The  assumption  of  rapid  relaxation  is  made 

12 
implicitly  in  the  analysis  of  Jackson  and  Nielsen.    If 

chemical  relaxation  occurs  too  slowly,  a  nonequilibrium 
layer  will  form  which  will  enhance  the  absorption  of  laser 
radiation. 

The  development  of  the  nonequilibrium  layer  occurs 
due  to  the  combined  absorption  of  thermal  radiation  and 
laser  radiation.  The  absorption  of  thermal  radiation 

induces  appreciable  ionization  ahead  of  the  spark  analagous 

31 
to  the  effect  of  precursor  radiation  in  shock  waves. -^   The 

free  electrons  are  then  selectively  heated  by  inverse 
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bremsstrahlung  absorption  of  laser  radiation  and  electron 
thermal  conduction.  The  electrons  quickly  relax  to  a 
Maxwell-Boltzmann  distribution  at  a  temperature  T^  as  a 
result  of  the  high  efficiency  of  electron-electron 
collisions  in  transferring  energy.  Since  the  transfer  of 
energy  to  the  heavy  particles  by  electron  collisions  is 
inefficient,  a  two-temperature  plasma,  with  electron 
temperature  exceeding  the  heavy-particle  temperature, 
develops.  Further  ionization  is  then  caused  by  collisions 
of  heavy  particles  with  high-energy  electrons.  Consequently, 
an  appreciable  electron  density  can  exist  in  a  layer 
ahead  of  the  spark  where  the  temperature  of  the  heavy 
particles  is  still  quite  low. 

The  inverse  bremsstrahlung  absorption  coefficient, 
given  by  Equation  (2.^)  in  Chapter  II,  is  proportional 
to  the  square  of  the  electron  number  density.  Therefore, 
the  development  of  a  nonequillbrium  layer  with  a  high 
electron  number  density  ahead  of  the  spark  means  that 
absorption  of  the  laser  radiation  and  consequent  heating 
-occurs  in  regions  where  the  heavy-particle  temperatures 
are  less  than  12,000  °K. 

Although  this  is  a  nonequillbrium  process,  the 
equilibrium  solution  obtained  in  Chapter  III  can  be  used 
to  illustrate  the  effects  of  the  anomalous  absorption  on 
spark  propagation.  Absorption  of  laser  radiation  in  regions 
with  temperatures  less  than  the  equilibrium  ignition 
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temperature  can  be  simulated  by  lowering  the  ignition 
thermal  potential  to  ©^  =  9-       <   .15  kW/cm.  The 
required  laser  intensity  for  a  given  propagation  velocity 
is  then  reduced,  according  to  Equation  (3.11),  to 


/ 


5.=      ^    S 

.15 


where  S  is  the  intensity  calculated  using  Q^    -   .15  kW/cm. 
Equations  (3.6)  -  (3.10)  indicate  that  the  spark  thermal 
potentials  are  reduced  by  the  same  ratio.  Therefore,  an 
absorption  wave  preceding  the  temperature  rise  lowers  the 
required  laser  intensity  and  the  maximum  spark  temperature 
for  a  given  propagation  velocity. 

For  example,  initiating  the  absorption  at  temperatures 
near  6,000  °K  (  9  c:   .025  kW/cm)  results  in  laser  inten- 
sities and  thermal  potentials  reduced  by  a  factor  of  six. 
Therefore,  a  free  spark  of  radius  .50  cm  propagating  at 
1  m/sec  would  require  only  123  kW/cm  and  would  result  in 
a  value  of  T^   =  17,000  °K.  For  an  ignition  temperature 
of  12,000  °K,  Sq  would  be  7^0  kW/cra^  and  T^^^  would  exceed 
4-0,000  °K. 

Figure  A-1  shows  the  isotherms  for  a  spark  of  radius 
.50  cm  propagating  at  1  m/sec,  assuming  Q^    =  .025  kW/cm. 
The  effects  of  convection  can  be  seen  in  the  blunt  spark 
front  and  elongated  isotherms.   It  is  interesting  to  note 
that  the  12,000  °K  isotherm  crosses  the  x-axis  at  x  =  .10  cm, 
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Since  the  6,000  \   isotherm  crosses  near  x  =  0,  the  width 
of  the  zone  in  which  the  temperature  falls  from  12,000  ^ 
to  6,000  °K  is  .10  cm.  This  distance  is  less  than  the 
absorption  length  of  .25  cm  for  the  vuv  radiation  emitted 
from  the  core.  Therefore,  the  model,  though  simple,  is 
at  least  consistent  in  that  the  radiation  can  penetrate  to 
regions  where  T  <  6,000  °K. 

The  intention  of  the  discussion  presented  above  was 
to  describe,  in  a  qualitative  way,  the  effect  re-absorption 
of  thermal  radiation  would  have  on  spark  propagation.  A 
future  investigation  keying  on  the  processes  occurring  in 
the  spark  front  would  be  useful  in  understanding  the  propa- 
gation mechanism.  In  particular,  a  one-dimensional  model 
of  a  nonequilibrium  gas  under  the  simultaneous  influence 
of  two  radiation  fields  would  be  a  logical  extension  to 
this  analysis. 
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